Parametric instability of linear oscillators with colored time-dependent noise 
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The goal of this paper is to discuss the hnk between the quantum phenomenon of Anderson local- 
ization on the one hand, and the parametric instability of classical linear oscillators with stochastic 
frequency on the other. We show that these two problems are closely related to each other. On the 
base of analytical and numerical results we predict under which conditions colored parametric noise 
suppresses the instability of linear oscillators. 

I. INTRODUCTION 

Let us consider the one-dimensional (ID) model described by the Schrodinger equation 

,-2 

-7^i^"{x) + U(xU(x) = EMx). (1) 
2m 

Here the ■)/'— function is a stationary solution for a particle of energy E moving in a random potential U{x). To simplify 
the form of the analytical expressions, in what follows we use energy units such that fi^ /2m = 1, and we set the zero 
of the energy scale so that the mean value of the disordered potential is zero, {U{x)) = 0. Here the angular brackets 
(...) denote the average over the disorder (i.e., over different realizations of U{x)). We restrict our considerations to 
the case of weak disorder, defined by the condition = ([/^(x)) <C 1. 

In the analysis of the model (1) one of the main questions is about global structure of the eigenstates ipix) in 
an infinite configuration space, — cxi < x < +oo. Of a particular interest is the problem of whether the eigenstates 
are localized or extended for x — > ±oo. As was shown in Ref. [1], in ID models any amount of disorder (even an 
infinitesimal one) results in the localization of all eigenstates (with the exception of a zero-measure set) provided that 
the potential U{x) is completely random. This means that the amplitude of every eigenstate decays exponentially 
towards infinity, therefore, far away from the localization center xq, one can write, 

\^ix)\ ^ exp{-\x-xo\/liE)). 

Here I (E) is the so-called localization length that characterizes, in average, the decrease of the amplitude of the 
eigenstate corresponding to the energy E. Analytical expression of I (E) that is relatively easy to obtain for a weak 
disorder, for which it is known that the localization length is inversely proportional to the square of the disorder 
strength, I ~ see below. 

Taking into account that the energy of a free electron is E ^ k'^ , the equation (1) can be written in the form of 
wave equation, 

i;"{x) + fcV(a;) = U{x)i^{x) (2) 

that describes wave propagation in different classical systems. One example is the propagation of electromagnetic 
waves in single- mode waveguides with a rough surface [2]. In this application the potential U{x) is determined by the 
horizontal profile ^(x) = e(p{x/Rc), where Rc is the correlation length of the profile, e <C d is the amplitude of the 
profile with d being the transverse size of the waveguide. Note that in this case the parameter k in Eq.(2) has the 
meaning of the longitudinal wave number k — ^ {'-^/c)'^ — (Tr/d)^ where u is the frequency of the wave. 

II. DISCRETE MODELS 

When the potential is constituted by a succession of delta scatterers, the model (2) takes the specific form 

n—oc 

r{x)^k'^{x)^ U,Mxn)5{x-x,,). (3) 
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Here Un is the amplitude of the nth delta-scatterer situated at x = x„. In experiments, potentials of this kind are quite 
easy to construct; in particular, one can obtain a realization of delta scattcrers by inserting an array of screws with 
predetermined lengths and appropriate positions in a single- mode waveguide [3]. Typically, one considers two limit 
cases. The first one occurs when all amplitudes C/„ are random variables, while the scatterers are periodically spaced, 
i.e., Xn = an . In this case one can speak of amplitude disorder. The second case is represented by the opposite 
situation in which the amplitudes of the scatterers are constant, [/„ — Uq while the positions a;„ are randomly 
distributed around their mean values, i.e., Xn = an + r]n with {rjn) = and (r?^) <C . Clearly, in the latter case 
the mean value of the potential is not zero; however one can handle this case within the framework of zero- mean 
potentials by making use of the special transformation to new variables, see details in Ref [4]. This second limit case 
can be referred to as positional disorder. 

Due to the delta-like form of a random potential, the model (3) can be considered as a discrete one. In fact, 
its analysis can be reduced to the study of an equivalent classical two-dimensional map which can be obtained by 
integrating Eq. (3) between two successive kicks of the scattering potential [4] , 

Pn+l = {Pn + An Qn) cos Un - gnSinfln, 

Qn+l = (Pn + An Qn) sin lln + 9„ COS /i„ . ' 

Here g„ and p„ are conjugate coordinates and momenta defined by the identities 

Qn = Tpn^ndpn = {^n COS/i„_l - 1pn-l) / Sin/X„_i 

where ipn is the value of the V'— function at the position x = Xn- The parameter /x„ is the phase shift of the ■^—function 
between two scatters, 

fin = k{Xn+i - Xn) (5) 

and the amplitude An of the nth kick is given by the the value of the potential at the position Xn, 

An = Un/k. (6) 

Free rotation in (4) between two successive kicks corresponds to free propagation between scatterers, and each kick is 
due to the scattering from a 6 spike of the potential. 

In the case of amplitude disorder, the phase shift between two successive scatterers is the same, /U„ = = ka, and 
the model (3) is known as the Kronig-Penney model. In this case the two-dimensional map (4) is equivalent to the 
following relation between -i/j^+i, ^„_i and ipn-, 

tpn+l + t/jn-l = COS 11+^ Sin/i^ Vn, (7) 

One can see that the relation (7) has the same form as discrete Schrodinger equation for the standard ID Anderson 
tight-binding model, 

^Jn+l + Ipn-l ^ {E + en) -ipn, (8) 

and describes electrons on a discrete lattice with the site energies e„. Therefore, many of the results for the Kronig- 
Penney model can be obtained by a formal comparison with the Anderson model, as discussed below. 



III. THE HAMILTONIAN MAP APPROACH 



One of the tools to find the localization length for discrete disordered models, cither analytically or numerically, is 
based on the transfer matrix method. In this approach the localization length can be expressed as the inverse of the 
Lyapunov exponent A which characterizes the growth of the eigenstates ^Pix) of the stationary Schrodinger equation 
for increasing x. An alternative approach can be obtained by interpreting the stationary Schrodinger equation as the 
equation of motion of a classical particle (in this scheme the space coordinate x of the disordered model is to be seen 
as the time coordinate for its dynamical counterpart). In particular, in the case of discrete disordered models, this 
approach leads to the study of classical maps. 

It is instructive to illustrate this approach by discussing its application to the simplest case of the Anderson model 
(8). Comparing Eq. (8) with Eq. (4), one can obtain that there is an exact correspondence between them by letting 
/U„ = ;U and 
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-E' = 2cos/u; j4„ = — e„/sin/U. 



(9) 



It is clear that for weak disorder the energy spectrum of the Anderson model (8) is close to the unperturbed one 
which is defined by the condition \E\ < 2; this legitimates the first equality in Eq. (9). 

To analyze the dynamics of the two-dimensional map (4), it is convenient to introduce the action-angle variables 
{rn,On) according to the standard transformation, q = rsinO, p = rcos6. As a result, the map gets the following 
form, 



where 



sin 6n+i = -D„ 1 (sin(6'„ - fi) - sin 6'„ sin fi) 

C0s6n+1 = D~'^ {cos{0n - IJ,) + An sin 9n COS jj,) 



Dn = ^ = ^l + A„sin(2^„)+A2sinX 



(10) 



Note that the following results for the localization length do not depend on the sign of n. It is important that the 
equation for the angle On can be written in the form of the one-dimensional map, 



C0t(^„+1 + fi) = cot 6n + An- 



(11) 



This fact simplifies the analysis of the distribution of On- The localization length I is defined as the inverse Lyapunov 
exponent, and the latter is determined by the standard relation [5] 
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Here the over bar stays for time average and the brackets for the average over different disorder realizations. The 
expression for can be splitted in two terms 



ri = (ln( 



> + (ln 



sin Qn 



+1 



sin On 



(12) 



The second term on the r.h.s. is negligible because it is the average of a bounded quantity. It becomes important 
only when the first term is also small, i.e. at the band edge |£^| « 2 or « 0. Thus, apart from this specific case, 
the localization length can be evaluated from the map (10) using only the dependence of the radius r„ on discrete 
time i„ = n. It is important to note that the ratio Vn+i/rn depends only on the angle 0„ and not on the radius r„; 
as a consequence, one can compute the localization length just by averaging the first term in (12) over the invariant 
measure p{9) associated with the ID angular map (11). 
In a direct analytical evaluation of (12) one can write, 

l-^ = j P{e) ln{D{e,e))p{e)d0de , (13) 

where P{e) is the density of the distribution of e„, and p{9) is the invariant measure for the angle variable. We 
use here the fact that p{9) does not depend on the specific sequence Cn, but can depend on the moments of P{e), 
particularly on its second moment tr^ (see discussion in [6]). As one can see, in order to evaluate the expression (13), 
first one has to determine the invariant measure p{0). 

In the case of weak disorder and not close to the band edges we have <C 1 and one can use the standard 
perturbation theory. This allows one to cast Eq. (13) in the form 
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2sin2yu 



2tt 

j e'P{e) de j p{9) Q - ^ cos(2^) + ^ cos(4e)) d9 . 



(14) 



This expression is valid for all energies within the band, but fails at the band edges, where one must include the 
contribution of the second term of Eq. (12) in the computation of the inverse localization length (see Ref. [6]). One 
can also see from Eq. (11) that in the limit of weak disorder the map for 9n has the approximate form. 



9n+l = On - P^ - An sin On + A^ sin ^On COS 



{mod 27r}. 



(15) 
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Therefore, in the first approximation the invariant measure for is flat and this makes possible an expHcit evaluation 
of the expression (14). Thus, one easily obtains 

= = (16) 

8 sin /X 8^1-^ 

where = (e^) is the variance of the disorder. It is interesting to note that the expression (16) is not correct at the 
band center, i.e., for = (see discussion and references in Ref [6]). The reason is that for this energy the standard 
perturbation theory fails and one must use specific methods to obtain the correct expression of As was found, 
the anomaly at the band center originates from the fact that for £ = the density p{9) is not flat, instead, it has a 
slight modulation with 6. This additional ^-dependence of the invariant measure is due to the special circumstance 
that the case E = corresponds to /i = tt/2 so that the map (15) becomes approximatively periodic of period four. 
As a consequence, p{6) has a weak modulation of period 7r/2 and, therefore, the fourth harmonic cos(40) in the 
expression (14) also gives a contribution. 

Due to the analogy between the Anderson model (8) (A-model) and the Kronig-Penney model (3) (KP-model) with 
fj,n = A*) 0116 can derive from the result (16) the expression for the localization length of the KP-model 



-\E) 



£q sin^(fca) 



8fc2 sin^ 7 

Here the phase 7 ( < 7 < tt) is given by the equation, 

2cos(fca) + — sin(A;a) = 2cos7. 
k 

This equation is the well-known dispersion relation for the periodic Kronig-Penney model; the parameter 7 plays the 
role of the Bloch number. 



IV. ANDERSON LOCALIZATION AND PARAMETRIC INSTABILITY 

It is easy to see that the Schrodinger eqiiation (3) for the quantum ID disordered model can be interpreted as 
the dynamical equation of a linear classical oscillator with a parametric perturbation constituted by a succession 
of delta- kicks at times tn = Xn- In particular, the map (4) corresponding to the Kronig-Penney model (/x„ = n) 
can be obtained by integrating the dynamical equations between two successive kicks for a stochastic oscillator with 
Hamiltonian of the form 

(22x2/^ \ 
y + yj + y E A^s{t-nT)\. (17) 
^ \n— — 00 / 

Therefore, g„ and p„ in Eq. (4) stand for the position and momentum of the oscillator immediately before the nth 
kick of amplitude An occurring at the time t = nT. Correspondingly, the phase shift between two successive kicks is 
given by /U = ujT where lo is the unperturbed frequency of the oscillator and T is the period between the kicks. 

In this description the exponential localization of the eigenstates of Eq. (3) corresponds to a parametric instability 
of the stochastic oscillator (17). The instability manifests itself as an exponential divergence of initially nearby orbits 
(orbit instability) and, correspondingly, as an exponential growth of the average energy of the parametric oscillator 
(energy instability). The Lyapunov exponent A, which gives the inverse localization length in the solid-state the 
divergence of classical trajectories (or, the rate of the energy growth). 

In the previous section we have considered the case of a weak uncorrelated disorder which is characterized by 
its variance only. In application to classical oscillators this corresponds to a white noise perturbation. In the 
following, we consider the general case of colored noise and show that noise correlations can lead to a quite unexpected 
phenomenon. To discuss the effects of correlated noise in parametric oscillators, we apply the approach of Ref. [7] to 
the continuous model described by the Hamiltonian, 

where ^(t) is a continuous and stationary noise. This model is slightly different from the one defined by Eq. (17) 
because the noise ^{t) is a continuous function of time rather than a succession of i5-kicks. We assume that the noise 
^(t) has zero average and that its binary correlator is a known function. 
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m=0 and m^{t + r))=x{r). 



(19) 



Here and below, in contrast to previous sections the symbol (. . .) will refer to the time average, {f{t)) = 

limTo^oo ^ Jq" f(t)dt, which is assumed to coincide with the ensemble average for the process ^{t). 
We define the Lyapunov exponent as follows, 

A= lim limlirin^.^. (20) 

As in the previous section, we introduce polar coordinates via the standard relations q = rsinO, p = rcosO. This 
allows us to represent Eq. (20) in the form 



1 Z"^" r 
A = lim — / - dt. 
To^oo To Jq r 



'0 

To proceed further, we consider the equations for the random oscillator in polar coordinates 

9 = Lo + ^{t) sin^ 9, 



r = -^r^{t)sm29. 

Using the last equation, the expression for the Lyapunov exponent can be finally written in the form 

A = lim r"msm{29{t)) = ^(^(t) sin (20(i))). (21) 

Therefore, the problem of computing the Lyapunov exponent (20) is reduced to that of calculating the noise-angle 
correlator that appears in Eq. (21). This was done in Ref. [7] by extending the procedure, originally introduced in 
Ref. [8] for discrete models, to the continuum case. As a result, the expression for the Lyapunov exponent takes the 
simple but non-trivial form, 

A = - / (C(i)C(i + r)) cos(2wr) dr. (22) 

One can see that the Lyapunov exponent for the stochastic oscillator (18) is proportional to the Fourier transform 
x{2co) of the correlation function at twice the frequency of the unperturbed oscillator. 

A similar result can be obtained for the parametric oscillator (17) with discrete noise. In this case the inverse 
localization length can be written as [8] 

A=-^<^(a;T). (23) 
Here the function (p {ojT) is the Fourier transform, 



+ 00 



(fi {uT) = 1 + 2 ^ C(A;) cos (2a;Tfc) 



fe=i 

of the binary correlator 

C(fc) = ^^f^ (24) 

of the colored noise. Therefore, the final expression is given by the product of two factors, namely, the Lyapunov 
exponent for the white noise case and the function (p(ujT), which describes the effect of the noise correlations (the 
color). In the case of white noise we have <f{u)T) = 1. 



5 



V. SUPPRESSION OF THE PARAMETRIC INSTABILITY 



Expressions (22) and (23) of the Lyapunov exponent for stochastic oscillators with weak frequency noise give a 
remarkable result: within the limits of the second-order approximation the rate of parametric instability depends 
only on the binary correlator of the noise. In application to solid state models this fact has suggested a way to 
construct random potentials with specific spatial correlations that result in "windows of transparency" in the energy 
spectrum. Indeed, if the Lyapunov exponent vanishes within some range of the energy (or, the wave number fc), then 
the corresponding eigenstates are extended in that energy interval. When one considers finite samples, this means 
that the transmission coefficient has to be one in the energy windows where the Lyapunov exponent vanishes. 

The possibility of engineering random potentials in order to obtain Lyapunov exponents with predefined energy 
dependence can be deduced from the expression (22) for continuous model, or from Eq. (23) for the discrete one. Both 
expressions show that if the Lyapunov exponent is known, the two-point correlator of the corresponding noise can be 
computed with an inverse Fourier transform. Since a stochastic process is not completely determined by its two-point 
correlator, one can conclude that there is actually an infinite set of noises which give rise to the same Lyapunov 
exponent because they have an identical binary correlator. 

As one can see, in order to have suppression of the parametric instability in classical oscillators with colored noise, 
one needs to have X{u)) = in some range of u). Although at first sight the construction of a random potential An or 
^(t) with a given binary correlator seems a difficult task, a rather simple method to solve this problem was presented 
in Ref. [8] for discrete models. This method was subsequently extended to oscillators with continuous noise. Here we 
describe how the method works for both classes of oscillators. Let us consider the continuous model (18) first. The 
starting point is the correlation function x(r) that can be easily obtained by inverting formula (22), 

Q l-OO 

x(t) = - / A(a;)e2-- ckj. 

Once the correlation function x(t) is known, we can obtain a stochastic process ^(t) satisfying the conditions (19) by 
means of the convolution product, 

/ + 00 
(3{s)ii{s + t)ds, (25) 
-oo 

where the function (}{t) is related to the Fourier transform x{u>) of the noise correlation function through the formula 

with r]{t) as any stochastic process such that 

(?7(i)) = and {'n{t)r]{t')) = S{t - t'). (26) 

Formula (25) defines the family of noises corresponding to a specific form X{u)) of the frequency-dependent Lyapunov 
exponent and gives the solution of the "inverse problem" (i.e., determination of a noise ^{t) that generates a pre-defined 
Lyapunov exponent). 

As an example, we consider the Lyapunov exponent 

^i^) ~ I ' otherwise ' ^'^'^^ 

whose frequency dependence implies that the random oscillator undergoes a sharp transition for = 1/2, passing 
from an energetically stable condition to an unstable one. Following the described procedure, it is easy to see that 
the Lyapunov exponent (27) is generated by a noise of the form 

X f^'^ sin(s) , , , 

^ i) = — / -^V{s + t) ds, 

with r]{t) being any random process with the statistical properties (26). Working along these lines, it is easy to see 
that one can also construct the frequency noises such that the parametric oscillator is stable for all values of u except 
those contained in a narrow frequency window. 
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We now turn our attention to discrete models of the form (17). To show how the expression that is equivalent to 
Eq. (25) can be worked out for this class of systems, we refer to the case of the Anderson model (8) with correlated 
disorder. Since this model can be put into one-to-one correspondence with the kicked oscillator (17), it is perfectly 
legitimate to analyse each of the two models in terms of the other; this approach has also the advantage of enhancing 
the physical understanding of the problem because it allows one to interpret the parametric instability of a stochastic 
oscillator in terms of localization of electronic states for the Anderson model. 

When we transpose the result (23) to the case of the Anderson model (8), we obtain that the expression for the 
localization length has the form, 

= (28) 
8sm n 

where 

oo 

^(^) = 1 + 2^C{k)cos{2ij,k). 

k=l 

Here ({k) is the binary correlator (24) which can be written in terms of the site energies e„ of the Anderson model as 

C(fc) = (e„6„+fe)/(e2). 

If the Lyapunov exponent (28) (therefore, the function (p{fi)) is known, the binary correlators (24) can be derived 
with an inverse Fourier transform, 

2 f^/^ 

C(fc) = - / ipiiJ,) cos(2/zfc) dii. (29) 
I" Jo 

As for the continuous model, the sequence of site energies e„ with the correlator of the specific form (29) can then be 
constructed with the convolution product. 



^ J2 f^kZn+k, (30) 



fe=— oo 



where 

2 r'^ 



2 ri^ 

Pk = - V'fiitJ') cos(2/iA:) dfj, 



and Zn are random numbers with the zero mean and unitary variance. It is easy to check that the correlators of the 
site potential (30) coincide with the Fourier coefficients (29). 

As an illustration of the method, we construct the random potential which results in the following function y(/i) 
for the Lyapunov exponent of the discrete Anderson model, 

_ / C'o if Ail < M < M2, 
^^^^ \ if < At < Ml or /Z2 < < 7r/2. 

Here, Cq = tt/2{^2 — A*i) is the normalization constant that results from the condition Co = 1- The corresponding 
localization length exhibits two sharp mobility edges at the values Ei = 2 cos a<i and E2 — 2 cos At2- Specifically, in the 
energy window Ex < E < E2 the eigenstates are strongly localized, while they are extended outside of this window. 
The binary correlators Q{k) for a random potential resulting in such a situation, are given by 

C(fc) = [sm{2k^i2) - suY{2kni)] . 
irk 

As a result, the expression for the inverse localization length reads 

2 

x = i- 



16 sin^ Ai [acos (-^) — acos (■^)] 
where cr^ = (e^). If the energy window A^; = Ex — E2 is narrow, one can write. 
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arccos 




One can see that the narrower the window A^;, the sharper the transition which occurs at the mobihty edges. This 
effect can be easily observed numerically, and it may have interesting applications for parametric oscillators. Indeed, 
small localization lengths correspond to large values of the Lyapunov exponent. Therefore, for values of the frequency 
of the kicked oscillator (17) which correspond to energy values within the localization window in the related Anderson 
model, the instability of the oscillator is very strong and one can speak of a kind of "parametric stochastic resonance" . 



VI. DISCUSSION 



In the previous section we have shown how a proper choice of colored noise (random potentials in the solid state 
models) can suppress the parametric instability of a stochastic oscillator in a prescribed frequency range. One should 
note, however, that the theoretical analysis has been focused on the case of weak noise and that almost all analytical 
results have been obtained using perturbative approach. Therefore, the conclusion that the Lyapunov exponent can 
vanish within some frequency region is valid only within the framework of second-order perturbation theory. Going 
beyond the second-order approximation, it is possible to estimate the correction to the present results for the inverse 
localization lenght and to show that the correction is represented by a term of order 0((T*) (with cr^ = (e^)). It is not 
clear whether one can make this fourth-order correction vanish with an appropriate choice of the statistical properties 
of the noise [9]; from a practical point of view, however, in the case of weak noise there is a well-defined separation 
between the time scale t ^ l/c^ over which the suppression of instability holds, and the much longer time scale 
t ~ 1/cr'' over which the effects of fourth-order corrections become relevant. When the second-order results for the 
inverse localization length are applied to wave-guides or solid state models, fourt-order effects can also be generally 
avoided by a proper choice of size of an experimental device [3] . 

The transfer matrix method in the Hamiltonian form described above is also very useful for finite times. In 
application to solid state models this question refers to transport properties through finite samples of size L. As is 
known [5,10], all transport properties can be directly related to the classical trajectories of the Hamiltonian map (4). 
Specifically, by studying general properties of these trajectories, one can find statistical properties of the transmission 
coefficient or the resistance. The transmission coefficient through a L-site sample can be expressed in terms of 
dynamical variables of the classical map (4) as 

J-T, — 7^ ^ 



where ri and r2 represent the radii at the Lth step of the map trajectories starting from the phase-space points 
Pi = (a^o = IjPo = 0) and = {xq = 0,po = 1), respectively. As for the resistance Rl, it is defined as the inverse of 
the transmission coefficient Rl = T£^. The key feature of these formulae is that they express the transport properties 
of a disordered sample in terms of the radii of map trajectories in the phase space. On the other hand, the square 
radius of a map trajectory is a quadratic function of the coordinate and momentum of the corresponding kicked 
oscillator, = p'^ + q^, and is therefore proportional to the energy of the latter. This fact makes possible to relate 
transport properties of quantum models with the time dependence of the energy of classical parametric oscillators. 

It is possible to obtain quite easily the moments of the energy of the parametric oscillator described by the 
Hamiltonian (17), see details in Refs. [7,11]. In particular, one can obtain that in the asymptotic limit (i.e., for times 
f 3> A^^) the mean value of the energy grows exponentially as 

{r^{t)) = r2(0)cxp(4At) 

where A is the Lyapunov exponent (23). This formula shows that the exponential rate of the energy growth for the 
parametric oscillator is four times the Lyapunov exponent, i.e., the rate of exponential separation of nearby orbits [7]. 

Another important question concerns the fluctuations of (t) for fixed times t depending on different realizations 
of the noise. Using the results of Ref. [7,11], one should distinguish between two different situations. The first one 
corresponds to small times when the value of r{t) is close to the initial value r(0). In solid state models this case is 
known as the ballistic transport for which the localization length I ^ is much larger than the size L = t oi the 
sample. At <C 1. Another limit case corresponds to large times, Xt^ 1, or to the strongly localized regime in quantum 
models. One of the most interesting effects is that in this case the fiuctuations of the energy of the classical oscillator 
(resistance in quantum models) are huge and the quantity i? = is not self-averaging. To deal with a well-behaved 
(that is, self-averaging) statistical property, one has to consider the logarithm of the oscillator energy, which has a 
Gaussian distribution for large times. It turns out that the energy has log-normal distribution: 
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P{r^t) 



exp 



(lnr2-2Ai) 



21 



8Xt 



This distribution implies that the energy of parametric oscillator, or the resistance R = of disordered samples, 
satisfy the relations 



(Ini?) = 2A; (In^ i?) = 4A + 4A'^ , 



Var(lni?) = (In^ i?) - (Inii)^ = 2{lnR). 

where A = Xt. 

In conclusion, we have discussed the analogy between properties of quantum ID models with random potentials and 
classical linear oscillators governed by parametric noise. We have shown that many results known for quantum models 
can be mapped unto corresponding properties of classical oscillators. One of the important questions is about the 
time-dependence of the energy of stochastic oscillators with frequency perturbed by a white noise. Another, even more 
exciting problem, is the behavior of the oscillators when the frequency noise has long-range correlations. It was shown 
that in the case of weak noise all statistical properties of the classical trajectories depend on the binary correlator 
of the noise only. This fact opens the door to the construction of colored noises with specific long-range correlations 
which result in a sharp change in the dynamical behavior of the parametric oscillator at some threshold value of 
the unperturbed frequency. Specifically, the characteristic instability of parametric oscillators can be suppressed in 
a certain frequency range (with a brisk transition), thanks to long-range temporal correlations of the noise. These 
results may find diff'erent applications in the field of classical systems with colored noise. 
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